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ABSTRACT: 

The model of heavy Wigner matrices generahzes the classical ensemble of 
Wigner matrices: the sub-diagonal entries are independent, identically dis- 
tributed along to and out of the diagonal, and the moments its entries are of 
order where TV is the size of the matrices. Adjacency matrices of Erdos- 
_ Renyi sparse graphs and matrices with properly truncated heavy tailed entries 

are examples of heavy Wigner matrices. We consider a family X^v of indepen- 
dent heavy Wigner matrices and a family Y^v of arbitrary random matrices, 
^ independent of Xjy, with a technical condition (e.g. the matrices of Yjv are 

I deterministic and uniformly bounded in operator norm, or are deterministic di- 

agonal). We characterize the possible limiting joint * -distributions of (Xjv, Yjv) 
in the sense of free probability. We find that they depend on more than the 
Ph * -distribution of Yjy. We use the notion of distributions of traffics and their 

free product to quantify the information needed on Yjy and to infer the lim- 
f~i iting distribution of (X^r, Yjv). We give an explicit combinatorial formula for 

^ joint moments of heavy Wigner and independent random matrices. When the 

d matrices of Y^r are diagonal, we give recursion formulas for these moments. 

I I We deduce a new characterization of the limiting eigenvalues distribution of a 

single heavy Wigner. 

I> keywords: *-distribution, asymptotic freeness, Wigner, heavy-tailed random variables, 

Erdos-Renyi graphs. 
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1 Introduction 
1.1 Motivations 

Let Ajv — {Ai, . . . , Ap) be a family of random N hy N matrices with complex entries, whose 
entries have all their moments. Following random matrix and free probability terminology, we call 
(mean) *-distribution of Aat the map 



r 1 

N 



Tr[P(A^)] 



defined on the set of non commutative * -polynomials, i.e. finite complex linear combinations of 
words in indeterminates ai, . . . , Op, aj^, . . . , a*. When it exists, we call limiting *-distribution of A^r 
the pointwise limit of $Ajv when N goes to infinity, and say that A^r converges in * -distribution. 



The notion of asymptotic *-freeness introduced by Voiculescu gives a rule to compute the 
limiting * -distribution of a large class of random matrices as their size goes to infinity (see 
m [nj [21 [T] [201 [m m] for examples). Recall its definition. 

Definition 1.1 (Asymptotic *-freeness). 

Let Ai, . . . , Ap be families of N by N random matrices having a mean limiting joint * -distribution 

P(Ai,...,Ap)] , 

defined on the set of non commutative * -polynomials in indeterminates ai,...,ap. The fami- 
lies Ai, . . . , Ap are asymptotically * -free if and only if for all * -polynomials Pi, P2, . . . , one has 
$(Pj(aij)) = 0, ij 7^ ij+i for all j ^ I implies (i>(Pi(ajJ . . . P„(ai^)) = for all n ^ 1. 

One of the main examples concerns independent Wigner and arbitrary random matrices. 

Recall that is a Wigner matrix whenever it is Hermitian with independent and centered sub- 
diagonal entries, such that the diagonal and the extra diagonal entries of yfNXjq are identically 
distributed according to probability measures, say v and /i respectively, that possess all their mo- 
ments. Let Xjv be a family of iV by independent Wigner matrices and Yjv a family of TV by 
N arbitrary matrices, possibly random but independent of Xjv. Assume that Yjv converges in 



Pv-^ lim E 

JV->-oo 



N 



;Tr 
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* -distribution and assume some control on Y^r (namely a concentration and a tightness property, 
see Assumptions 2 and 3 below, e.g. the matrices of Y^r are deterministic and uniformly bounded 
in operator norm). Then Voiculescu's asymptotic freeness theorem [1, states that the limiting 
* -distribution of (Xat, Yat) exists and depends only on the limiting * -distribution of Ypf and of 



the variances of extra diagonal entries of the Wigner matrices (see Corollary 3.9 1 



This fact reflects a universality phenomenon for eigenvalues statistics of large random matrices, 
since the limiting *-distribution of (Xat, Yn) does not depend on the details of the law of 
the entries of Wigner matrices. Such a result is useful since the convergence in * -distribution 

of (Xjv,Yjv) implies (and is actually equivalent to) the convergence in moments of the (mean) 
empirical eigenvalues distribution of any Hermitian matrix Hpj, obtained as a fixed *-polynomial 
in Xjv and Y^y. Recall that the empirical eigenvalues distribution of an N by N matrix Hpj is the 
probability measure 



1 ^ 



where Ai, . . . , \n are the eigenvalues of Hj^ and 5\ denotes the Dirac mass in A. 

This theorem of asymptotic *-freeness has an impact in classical probability theory. Consider 
the case where the families Xjv and Y^r consist only in one Hermitian matrix Xj^ and respec- 
tively. Then, the convergence in * -distribution of is the convergence in moments of its empirical 
eigenvalue distribution toward a probability measure tt. Assume the technical conditions stated in 
Assumptions 2 and 3, Section |3.1[ Then, Voiculescu's asymptotic freeness theorem gives in par- 
ticular a characterization of the limiting empirical eigenvalues distribution of the sum Yj^ -I- X^y. 
Following free probability terminology, it is called the free convolution of tt with the semicircle 
distribution a a, namely the probability measure 

daa{t) = 77^^ - r = 2Va= lim 2(E[7VXjv(1, 2)^]) ^ 

In this article we investigate the question of the convergence of (Xjv, Y^r) where Wigner matrices 
are replaced by matrices of a larger class, with the same structure of independence of its entries. 

Definition 1.2 (Heavy Wigner matrices). 

A random matrix Xff is an N by N heavy Wigner matrix whenever: 

1. Almost surely, Xiq is Hermitian, i.e. X^^ — 

2. the sub-diagonal entries of Xn are independent and centered, 

3. the diagonal entries of — y/NX]\j are distributed according to a measure vjq on M, 

4.. the strictly sub-diagonal entries of Mn are distributed according to a measure fj.^ on C, 
invariant by complex conjugacy. 



5. fiN and i^n possess all their moments and for any 1 



J t^'^diypfit) 



ak, (1.1) 
0(1). (1.2) 



The sequence {ak)k^i is called the parameter of X^. Hence, a Wigner matrix is a heavy Wigner 
matrix such that, with the notation above, the measures vn and /ijy do not depend on N. In that 
case, Ofc is zero for any fc ^ 2. Such a parameter is said to be trivial in the following. 



This matrix model has been introduced independently by two authors. Zakharevich [25] has 
shown that the empirical eigenvalues distribution of a heavy Wigner matrix converges as N goes 
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to infinity. She has proved that the Hniiting distribution depends only on the parameter of the 
matrix. It consists in the semicircular distribution of radius y/oi if the parameter is trivial, and has 
unbounded support otherwise. Zakharevich has given a combinatorial formula for the moments of 
this limiting distribution, based on the enumeration of certain rooted trees and she has proved that 
these moments characterize the measure when ak — 0{a'^) for some a > 0. Furthermore, Ryan \22\ 
has proved that a family Xjv of independent heavy Wigner matrices has a limiting *-distribution. 
He has given a combinatorial formula for it, which involves partition generalizing the classical ap- 
proach for large Wigner matrices based on non crossing pair partitions (see [21 ). In particular, he 
has proved that Voiculescu's rule of *-freeness does not govern the limiting * -distribution of Xat 
as soon as at least two matrices of the family have a non trivial parameter. Motivated by question 
from free probability, Benaych-Georges and Cabanal Duvillard [5 have shown the convergence 
of the empirical eigenvalues distribution for the generalized Gram matrix Hj^ m — Xn,M'^*n 
where Xn^m is an N by M matrix such that: 

• the ratio ^ converges to a positive constant 



the entires of v MXn,m are independent and identically distributed according to a probability 



measure /i^v on C whose moments satisfies Ryan-Zakharevich's condition (1.1 1 



In this article, we consider a family Xjv of iV by independent heavy Wigner matrices and a 
family of N by N matrices, possibly random but independent of Xjy. We characterize the 
possible limiting * -distribution of (Xjv, Yjv) under suitable assumptions on Y^. The most mean- 
ingful phenomenon that arises is that the limiting *-distribution of (Xat, Yat) depends on 
much more than the limiting *-distribution of Yjy. We use the notions of distributions of 
traffics and their free product to specify asymptotic statistics on Y^r and then characterizes the 
* -distribution of (X^r, Y^r). 

In particular, if Xj^ is a heavy Wigner matrix and a random Hermitian matrix, independent of 
Xj^ having a limiting empirical eigenvalues distribution and uniformly bounded in operator norm, 
the problem of characterizing the limiting eigenvalues distribution of Y/v +Xpj is ill-posed. Up to a 
subsequence, a limiting eigenvalues distribution exists but there can exist many possible limits. For 
instance, if Yn is a Wigner matrix independent of Xn, then the limiting eigenvalues distribution of 
Yn + Xn is the free convolution of Zakharevich's distribution with a semicircular distribution. If 
Y]^ is an arbitrary matrix with fixed limiting eigenvalues distribution (say diagonal and random), 
then the one oi X^ + Y^ depends on the whole distribution of traffics of Y^y- 

1.2 Examples of models 

We point out how the study of heavy Wigner matrices is interesting since this model is related to 
classical random matrices. 

1.2.1 Matrices with truncated heavy tailed entries 

We say that a law of a random variable x belongs to the domain of attraction of an a stable law 
if there exists a function L : M — > M slowly varying such that 

P(|x| ^u)^ ^^,Vw e M, a e]0,2[. 
\ / yd 

A Levy matrix Xpj with parameter a in ]0, 2[ is a random Hermitian matrix such that: for any 
i,j = l,...,iV, 

where the random variables (a;i,j)i^i^j^Ar are independent, identically distributed according to a 
law that belongs to the domain of attraction of an a stable law for an a in ]0, 2[ and 

t7Ar = inf{ ueR+\ P{\xi,i\^u) ^ ^ }• 
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By the formula for truncated moments of heavy-tailed random variables |31 Formula (15)], for any 
B > 0, the random matrix whose entries are given by: for any i, j = 1, . . . , N, 



-N 

is a heavy Wigner matrix with parameter 



" ^ (1.3) 



V(2/c - a)B°'Jk^i' 

The first mathematical results on Levy matrices are due to Ben Arous and Guionnet [4 in 2007, 
who have shown the convergence of the eigenvalues distribution of a single Levy matrix. Belinschi, 
Dembo et Guionnet [3J has studied the limiting spectrum of the sum of a Levy matrix and a 
diagonal matrix, and of a band Levy matrices. Moreover, Bordenave, Caputo and Chafai' [6^ has 
given an other characterization of the limiting distribution of a Levy matrix than one in [4J. It 
is based on the local operator convergence of a Levy matrix to a certain graph whose entries are 
labelled by random variables, the Poissonian weighted infinite tree. This is reminiscent with the 
traffic based approach of this paper for heavy Wigner matrices. 

1.2.2 Adjacency matrices of graphs and networks 

Let Gm = {y^E) be a simple undirected random graph with N vertices labelled {!,... ,-/V}. The 
adjacency matrix of Gat is the matrix 



An — (l{TO,n}e£;) 



.N' 



By simple, we mean without loops nor edges, so that is a symmetric matrix with entries in 
{0, 1} and its diagonal elements are zero. 

Erdos-Renyi sparse graphs: The only random graph, invariant by re-indexation of its ver- 
tices and whose adjacency matrix has independent entries is the Erdos-Renyi random graph: it is 
a random undirected graph with vertices {1, . . . , N}, such that two distinct vertices are linked by 
an edge with probability pM, independently of the others edges. 

We consider piq of the form for a fixed a > and N large, and denote cj^^ a random Erdos- 
Renyi with that parameter. This is called the sparse regime. 

Consider the adjacency matrix A^-* of and set 

^(c) _ Ac.) a 



where is the N hy N matrix with zero on the diagonal and one elsewhere (X^r is made in such a 

Ac.: 



way its entries are centered). Then, A^^"* is a heavy Wigner matrix with parameter {a)k^i. It can 



be observed that A^'' and A^'' have the same limiting eigenvalues distribution. More generally, 
one can replace heavy Wigner matrices by adjacency matrices of sparse Erdos-Renyi graphs in the 
results of this article. 

Network version: With A^'' as above, denote by A^-* the Hermitian random matrix obtained 

from j4^^ by replacing its non zero entries by independent, identically distributed random vari- 
ables. More formally 

X^^^^A^^KMm, 

where o denotes the Hadamard (entry-wise) product, and Mjv is a random Hermitian matrix, 
independent of Aj^' , with independent identically distributed entries (up to the Hermitian condi- 
tion) . Assume that the common law of the entries of are distributed according to a measure 
centered and which possesses all its moments. Then, A^^ is a heavy Wigner matrix with parameter 

(ax 1 |zp^-dM(^))^^^. (1.4) 
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The spectral theory of weighted graphs is a fields that has been intensively developed in mathe- 
matics [nilE]. The analysis of the spectrum of random sparse matrices has been tackled by many 
authors as Khorunzhy, Shcherbina and Vengerovsky [2] , Ding and Jiang [TU] and Shcherbina and 
Tirozzi 

Remark that, more generally, if Xpf is a heavy Wigner matrix with parameter {ak)k^i and Afjy is 
as above, independent of Xn, then o Mn is a heavy Wigner matrix with parameter 



(1.5) 



2 Statement of the main results 

2.1 The limiting distribution of independent heavy Wigner and arbi- 
trary random matrices 

Notations: Let Yjv — (Y^^\ . . . ,Yq^^) of N hy N random matrices. Let Un be a uniform per- 
mutation matrix, independent of Yjv and denote Zj — UnYj^^'U'^ for any j — 1, . . . ,p. Following 
Lovasz [15 , we call injective density of T in Y^v the quantity 



r K 



S%[T{Yn)] =E 



Y 



N 



(2.1) 



J = 1, 



is a set of indices 



where T = {(fc„Z,,7(j),e(j)) e {1,2,...}2 x {l,...,p} x {1,*} 

(fixed as N go to infinity), M{k, I) denotes the k, I entry of a matrix M, M* its conjugate trans- 
pose, and E[ • \Y n] means the conditional expectation with respect to Yat. 



Such a set of indices T is seen as a labelled graph: with the notations above, the set of ver- 
tices is \^ = {kj, Ij I j — 1, ... , K} and the multi-set of edges is i? = {{(^jj Ij) \ j — ■ ■ ■ ^ ^}}- 
Each edge e has a label If T is connected, we call it a *-test graph. It is called a cyclic 

*-graph when the integers (fcj, Ij) can be taken of the form (kj, fcj+i), j — 1, . . . , K with kx+i = ki. 

We state three assumptions on a family of random matrices. 

Assumption 1: One assumes its convergence in distribution of traffics on Gcyc{y,y*) (see 



Section 3.1|, that is, the convergence of the quantities 



E 



E 



{N-\v\y. 



for any cyclic * -graph T with \V\ vertices. This mode of convergence extends the convergence in 
* -distribution and the weak local convergence of graphs [IB]. 

Assumption 2: One assumes a concentration hypothesis in this setting: for cyclic *-test graphs 

Tl, ■ ■ ■ , Tn , 



E 



r^[r„(Y^)] 



E 



N 



..E 



In 



[TniY 



N 



0. 



Assumption 3: One assumes a condition which implies tightness in the setting of the main 



theorem of this paper (Theorem 2.1 Section 3.2.2): for any (non cyclic) *-test graphs Ti, . . . ,T„ 



E 



r^[Ti(Y^)] ...rj^[r„(Y^)]] =0(7V^"-W^-)/2-i 



where r(Ti) is a positive integer defined by Mingo and Speicher P^, called the number of leaves 
of its tree of two-edge connected components. It depends on the geometry of Ti, see Section fS.l. 2 [ 
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We can state the main theorem of this article, where the interest in that it gives an expHcit 
description of the hmiting objects (by the traffic freeness, whose definition is recalled in Section 



Theorem 2.1 (The traffic-asymptotic freeness of Xi, . . . , Xp, Yjv). 

Let be a family independent heavy Wigner matrices, independent of an arbitrary family Y^r 
of random matrices satisfying the three assumptions above. Then, the joint family (XtvjYjv) has 
a limiting * -distribution. Its limit is characterized by the notion of traffic-freeness in the sense of 
fl&j . and depends only on the parameters of the heavy Wigner matrices and of the distribution of 
traffics o/Yat on t/cyc(x, x*). 



Theorem |2.1 1 is proved in Section f3. 2. 2[ In the case where the matrices of X^r are classical Wigner 
matrices, we can replace Assumption 1 by the convergence in * -distribution of Y^r (see Corollary 
3.9 1. To the author's knowledge, this improves the usual asymptotic freeness theorem for indepen- 
dent Wigner and arbitrary random matrices. 



We give examples of families of matrices which satisfy assumptions of Theorem |2.1[ 
Proposition 2.2 (Examples of matrix models). 

1. Matrices uniformly bounded in operator: Let Y ^ be a family of random matrices 
whose operator norm is almost surely uniformly bounded. Then, up to a subsequence, Y^r 
satisfies Assumption 1 and it always satisfies Assumption 3. 

2. Diagonal matrices: Let Ti^ be a family of random diagonal matrices having a limiting 
* -distribution and satisfying: for any * -polynomials Pi, ... , P„, 



E 



1 



Tr[Pi(DAr)]...-Tr[P„(D^)] 



E 



Tr[Pi(D^)] 



.E 



Tr[P„(DM)] 



0. 



Then, D^v satisfies the three assumptions of Theorem\2.1 



3. Adjacency matrices of graphs: Let Yjv be the family of adjacency matrices of a random 
colored graph Gn, with uniformly bounded degree, that converges in the sense weak local 
convergence to a random rooted colored graph. Then, Y^r satisfies Assumption 1. Moreover, 
if it satisfies Assumption 2, then it satisfies Assumption 3. 

As we can use this theorem for Y^r containing diagonal matrices of projection, we obtain an 



analogue of Theorem 2.1 for covariance matrices (see Proposition 3.101. We also obtain the weak 



convergence of the empirical eigenvalues distribution of Hermitian matrices in independent Levy 



and random matrices (see Proposition 3.11) 



2.2 The limiting distribution of independent, permutation invariant ran- 
dom matrices 



Theorem |2.1| is a consequence of a more general theorem, namely the traffic-asymptotic freeness of 
random matrices on cyclic *-test graphs stated below, which specify a result of in the settings 
of the three assumptions stated above. It is proved in Section 3.2.2| once the setting of traffics and 



their free product have been reminded in Sections |3.1| and |3.2.1 



Theorem 2.3 (The asymptotic traffic-asymptotic freeness of permutation invariant, independent 
families of matrices) . 

Let . . . , Zj^'* be families of random N by N matrices. Assume that 

• Z^-*, . . . , Z^-* are independent, 

• for any j — 1, . . . ,p, Z^-* is permutation invariant, except possibly for one j in {1, . . . 

• Zj^' satisfies the three assumptions of Theorem [1.2]. 
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Then, the joint family (Z^^'' 



N ' ■ 



Z^-*) has a limiting distribution of traffics on Gcyc{z,z*) , which 



is the traffic-free product of the limiting distributions of Z 



(1) 

N ' ■ 



r(p) 
-'N 



on g. 



cyc\ 



In par- 



ticular, it has a limiting * -distribution which depends only on the limiting distribution of traffics 
of Zj^\ . . . , Zj^-* separately. 

We prove that a single lieavy Wigner matrix satisfies tliese assumptions. Hence, consider indepen- 
dent heavy Wigner matrices Xjv — {x[^\ . . . , Xp^'') and random matrices satisfying the three 
assumptions. We get that the families (xj^-* )...., (Xp^''), Yjv are asymptotic trafhc-free (see 
Section 3.2.1| as N goes to infinity. The rule of trafhc-freeness gives an explicit way to compute 
*-moments in (Xjv,Yjv). It is based on the computation of the limiting distribution of traffics 
given in the following proposition. 

Proposition 2.4 (The limiting distribution of traffics of a heavy Wigner matrix). 
Let Xn be a heavy Wigner matrix with parameter (aj,)fe^i. Then, Xjq has a mean limiting distri- 
bution of traffics on Q{x) given by: for any cyclic * -test graph T, 



E 



r%[T{Xr,)\ 



n 







ifT is a fat tree of type {qk)k^i, 
otherwise. 



A fat tree is a * -test graph which becomes a tree if we forget the multiplicity and the orientation of 
the edges. A fat tree is of type {qk)k>i if it has qk undirected edges of multiplicity 2k. See Figure 

m 



C 
Z 






Figure 1: Left: the first fat trees with two vertices. From left to right, they are of type 
(1, 0, 0, ... ), (0, 1, 0, 0, ... ) and (0, 0, 1, 0, 0, ... ). Other examples: from left to right, two fat trees 
of type (4, 0, 0, ... ) and one of type (2, 1, 1, 0, 0, . . . ). 



2.3 Limiting joint *-moments of (X7v,Y 



N) 



Based on the results of Section 2, we give a combinatorial description of the limiting * -distribution 
of (Xjv, Yjv). This approach is different that Ryan's one by the so-called clickable partitions. 



It can be considered as a dual version (see Figure 10 1. This generalizes the description |12j for 



non-heavy Wigner matrices by rooted, oriented trees and the description [T3] for a single adjacency 



matrix of a Erdos-Renyi sparse weighted graph (see Section 1.2.2) by minimal walks on such trees 



We first precise the vocabulary we use in order to avoid ambiguities. A tree is a graph with 
no cycles. We call rooted tree in the complex plane a undirected tree possessing one marked vertex 
(called its root) embedded in the non negative half plane of by planting its root at the origin. 
A directed edge of such a tree refers to a pair of adjacent vertices. A cycle on a tree is a sequence 
of directed edges of the form (^{vi,V2), . . . , {vl-i,vl), {vl,vi)). With this notation, L is called the 
length of the cycle and the directed edge (d„, Vn+i) is called the n-th step of c (n = 1, . . . , L with 
indices modulo L). 

From now, we fix a *-polynomial P — x~f{X)P\(y) . . . x^i^jj^P-^i^), where L ^ 1, 7 : {1, . . . , L} — >■ 
{!,... and Pi, . . . ,Pl are polynomials. We set 

a>(P)= hm E[i-Tr[P(X^,Y^)] 

N ^00 I\ 
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and give a combinatorial description of this quantity. By linearity and traciality, this characterizes 
the limiting * -distribution of (X^r, Yat). 

Definition 2.5 (Colored, minimal cycles on trees). 

We set C^^^ the set of all couples {G,c), where G is a rooted tree in the complex plane with less 
than [^J edges, and c is a cycle on G with the following properties: 

• c starts at the root of the tree. When it visits a new vertex, it visits the leftmost one. It visits 
all the vertices of G and has L steps. 

• By convention, we say that the n-th step of c is of color ^{n), n ^ 1, . . . , L. Then, c must 
visit each edge of T with a single color. 

We have drawn some examples of minimal cycles on trees (in one color) in Figure [2] From now we 
fix (G, c) in and define separately weights associated to heavy Wigner and other matrices. 






Figure 2: Examples of minimal cycles on trees (in one color). The root of the tree is marked with 
a black dotted point. The start of the cycle is announced by a black square. We figure out the 
direction of the cycle each time it approaches a vertex. 



Definition 2.6 (Heavy Wigner weights). 

For any edge e of the tree, we denote by j(e) the color in {1, . . . ,p} of steps on this edge, and by 
2k{e) the number of times the cycle visits it. We set 

l^HwiG.c)^ Y[ aj(e),fe(e)) (2-2) 
e edge of G 

where for any j — 1, . . . ,p, {aj^k)k^i stands for the parameter of Xj'^K 
Definition 2.7 (Traffic weights). 

Write the cycle c = (ei, . . . ,6^), where Cj is a directed edge of the tree. For any vertex v of the 
tree G, we define a * -test graph Ty in the variables Pi(y), . . . ,PL(y)- The vertices of Ty are the 
incident edges of G in v. If the n-th step of c is incident at v, then we get an edge between the 
undirected edges corresponding to e„ and e„+i (with the convention e^+i = ei) which is labelled 
F„(y). We set 

^™(G,c)= H T[Ty], (2.3) 

V vertex of G 

where for any * -test graph T, with set of vertices V , set of edges E and whose edge e in E is labelled 

e{e) 

4>:V^{l,...,N} eeE 



(see Section 3.1 ) 
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Figure 3: Left: A colored, min- 
imal cycles on a tree. Right: 
We extract the features of the 
cycle by drawing boxes around 
each vertex and forgetting the 
initial tree. The edges out of the 
boxes give the weight lohw'- we 
count the number and the color 
of visits of each edges and use 
formula (2.2 1. On the right of 



each box, we have represented 
the associated *-test graph in 
variables Pi(y), . . . , PL(y). The 
product of their trace of *-test 



graph (see Section 3.11 gives the 
weight lotr- 







Figure 4: Some boxes that can be obtained by the process of Figure [3] and the graphs obtained 
from them. 



See Figures [3] and |4] for illustrations and concrete procedure to compute these weights. 

Proposition 2.8 (Joint moments in heavy Wigner and deterministic matrices). 
For any polynomial P of the form P = a:^(i)Pi(y) . . . Xj(^]^^PL{y), one has 

$(P)= J2 ^HwiG,c)XLJTB.{G,c). (2.4) 
(G,c)e£(T) 

We deduce from Proposition |2 . 8 1 the following simple facts. 
Corollary 2.9 (Basic properties of limiting variables). 

1. For any j — 1, . . . ,p, ^{xj) ~ and $(a;|) = CLj,!- 

2. For any L ^ \, any 7 : {1, . . . , L} — > {1, . . . ,p} and any * -polynomials Pi(y), . . . , Pl{y)) the 
quantity <i>(a::-y(i)Pi(y) . . ■ a;^(i)Pi(y)) vanishes as soon as the number of occurrence of one 
variable is odd. 

3. for any integers ni,...,?!^ ^ 0, any distinct indices in {l,...,p} and any *- 
polynomial P, one has . . . a;"/'P(y)) = ^(x"^^) . . . $(x"/')$(P(y)) . 

2.4 Applications 

2.4.1 The non asymptotic *-freeness of heavy Wigner and independent arbitrary 
matrices 

The notions of *-freeness and traffics freeness are different (TB]. In particular, by Definition |1.1| if 
the matrices Xn and are asymptotically *-free, then we get by definition |1.1| 

/(x, y) := - <^{x^)) (y2 _ $(^2)) (^2 _ ^^^2)^ (y*2 _ <j,(y*2))^ ^ ^ 

The following Proposition extends the result of Ryan which states that independent heavy 
Wigner matrices with non trivial parameters are not asymptotically *-free. We define the bilinear 
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forr 



(^"1,^2) lim E 



-Tr[Pi(XAr,YAr)oP2(Xjv,Y 



N) 



(2.5) 



defined for non commutative *-polynomials Pi and P2, where o stands for tlie Hadamard (entry- 
wise) product. This quantity is well defined since it can be written as the limit of the trace of 
cyclic *-test graphs in X^v, Y^v when Pi and P2 are monic *-monomials (see |16)). 

Proposition 2.10 (The non asymptotic freeness of heavy Wigner and deterministic matrices). 
Let X]\j be a heavy Wigner matrix with parameters {ak)k^i. Let Y/v be an arbitrary random matrix, 
independent ofX^ and satisfyin g the assumptions of Theorem \2.1\ Denot e by $ their joint limiting 
* -distribution given by Theorem 2.1 and the bilinear form given by (2.5). Then, one has 

f{x,y)=a2 X g{y), g{y) = '^'^^\y'^ ^y*"^) - |$(y^)|^ 



Remark that 02 is nonzero as soon as the parameter of Xj^ is not trivial (Proposition A.ll. Here 
are examples of matrices Yat as in the theorem such that (X^r, Y^y) are not *-free. 

Proposition 2.11 (Example of random matrices non asymptotically *-free from Xat). 

L IfY^ is a heavy Wigner matrix with parameter {bk)k^i, then g{y) — 62. 

2. IfY^ is a diagonal matrix having a limiting * -distribution, then 



1 



1 



g{y) - lim E -Tr[Y^Y^'] - -Tr[Y^] 



N 



.N 



N' 



Hence, (Xjy, Yjy) are not asymptotically * -free as soon as the limiting eigenvalues distribution 
of Yn is not a Dirac mass. 



2.4.2 A system of Schwinger-Dyson equation for the limiting distribution of inde- 
pendent heavy Wigner and diagonal matrices 

We prove recurrent relations for the joint moments in (X^r, Y^f) in the case where the matrices 
of Yjv are diagonal and it satisfies Assumptions 1 and 2. The philosophy of the proof is the same 
as for the recurrence relation on sparse, weighted, random graphs |14| . A difference with our 
approach is that these equations involve only moments in the entries of the matrices rather than 
purely combinatorial quantities. 



For any integer K ^ 1, we set the X-linear forms 

1 



: {Pi,...,Pk)^ lim E 



N 



Tr [Pi (Xjv, Yjv) o . . . o Pk(Xn,Yn ) 



defined for *-polynomials Pi, ... , Pr- in variables x, y, where o denotes the Hadamard (entry-wise) 
product of matrices. The maps ($^^')/<-^i are simple examples of statistics of distributions of 
traffics that are not defined for * -distributions. 

Theorem 2.12 (A Schwinger-Dyson system of equations). 

For any j = l,...,p, we set {aj^k)k^i the parameter of the matrix X^^K Then, the family of 
linear forms satisfies the following equations. For any integer K ^ 1, any monomials 

Pi, ... , Px and any j = 1, . . . ,p, one has 



$(^)(x,Pi,P2,...,Pif) = ^a,,fc 

fc>i 



E 

\-SK—k L.R. 



^$W(L)$(fc+^-i)(R), (2.6) 



where the last sum is over all the families of monomials 



L = (4'\-.-,4l 



-,(1) 



R — (Pi^\ . . . , P^^\ Pq 



(2) 



P(2) 



, R. 



P(^)) 
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such that 



Pk = Ri''\x,L[''^Xj)R[''^ . . . {xjLi'l^x,)Ri1\ k = 2,...,K. 

This gives a characterization of the hmiting * -distribution of (X.pf,~Ypf) and a way to compute 
joint * -moments. 

2.4.3 The spectrum of a single heavy Wigner matrix 

We deduce from Schwinger-Dyson system of equations a characterization of the spectrum of a 
single heavy Wigner matrix X]\j. Denote by (afe)j,^i its parameter. For any _ftr ^ 1, we set the for 
formal power series G^{K) in ^ 

g'w-^Ef ^ <i>(^)(x-,...,^-). 

n^O Tii + ...nK=n 

This quantity is simply a formal analogue of 

$(^)((A-:rr\...,(A-x)-l), 

as the terms (A — x)^^ are expended into formal power series (A — x)^^ — ^ X]n>o 5^^-^"' ^'^^ ^^^^ 
order between and these sums are interchanged. In particular, G^(l) is a formal analogue 

for the limit of 

E 



lTr[(A-X^.)-i] 



called the Stieltjes transform of X 



N 



Proposition 2.13. 

For any K ^ 1, we have the equality between formal power series in j 

G\K) = l(^G\i^-l) + 5]a,.(^^ + ^^2)G^(fc)G^fc + i^-l)). (2.7) 

These equations characterize the sequence {G^{K)^^^^ among the set of formal power series 
{G^{K)^j^^^ such that for any K ^ 1, the valence of G^{K) is larger than K. 

Remind that if we assume that there exist some constants a,/? > such that for any k ^ 1, 
flfc ^ afc'^, then by a result of Zakharevich |25| . the limiting eigenvalues distribution of is 
characterized by its moments. 
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3 The convergence of heavy Wigner and arbitrary matrices 

We first remind some definitions and results in [16 . We recall and reformulate the two first 
assumptions of our main theorem. 

3.1 The distribution of traffics of matrices 
3.1.1 Statement of Assumptions 1 and 2 
Definition 3.1 (*-test graphs). 

1. A * -test graph in variables x = (xi, . . . , Xp) is a finite, connected, oriented graph (with 
possibly multiple edges and cycles) whose edges are labelled by variables xi, . . . , Xp, x\, . . . ,x*p. 
Formally, it consists in a triplet T = [V, E, 7, e) where (V, E) is a graph and 7, e are maps 
from E to {1, . . . ,p} and {1, *} respectively, in such a way an edge e in E is labelled x'l!^(^^y 

2. A* -test graph is said to be cyclic whenever there exists a cycle on its graph visiting each edge 
once in the sense of its orientation. We denote Qcyc{'x.,x*) the sets of cyclic * -test graphs in 
indeterminates li. (we keep these notation even though the notation for indeterminates can 
change). 

Let Yjv = (^1, • ■ • ,Yp) be a family of by random matrices. For any *-test graph T. we call 
the trace of T in Y^v the quantity 

Traffic moments: TAT [T(Y^)] := ^ ^ R^tSW^))' (^-^^ 

<p:V^{l,...,N} eeE 

where 

• for any directed edge e = {vi,V2), we have set 0(e) = (0(ui), 0(w2)), 

• and for any N hy N matrix M and any integers n, m in {1, . . . , N}, the complex number 
M{n, m) is the entry (n, m) of M. 

• M* is the conjugate transpose of the matrix M. 

With the same notations, we call the injective trace of T in Y^v the quantity 

Traffic cumulants: TO[r(Y^)] := 1 ^ H ^^t) ('^(^))- (^.2) 

4>:V^{l,...,N} e£E 
injective 

The definition of is consistent with the definition of the introduction. Indeed, let gn be a 
random permutation of {1, . . . , N} independent of Yjy and JJ^y the random permutation matrix 
associated to ctjv- Then, for any *-test graph T, 



Y 



0:V-^{1,...,7V} ^eeE 
injective 

I E A\{{u.Y^Uu*,){m 



N 



N 

<t>:V^{l,...,N} 
injective 



Y 



N 



'"-''VMT(Y«,]. 



{N-\V\)< 



where ^^v the injective density defined in (2.1 1 



One can compute the trace of *-test graphs in terms of injective traces and vice versa. This 



is a consequence of Formula (3.3 1 stated below and of simple facts on posets (see HU). 
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Proposition/Definition 3.2 (Trace and injective trace of *-test graphs). 

Let T be a * -test graph whose set of vertices is denoted by V. Let tt be a partition ofV. We denote 
by T^ the * -test graph obtained by identification of vertices that belong to a same block ofir, see 
Figure [5| 

1. Matrix setting: For any family Ypf of N x N matrices and any * -test graph T, 

t4t{Y^)]= J2 A[T^Ym)], (3.3) 

where V{V) is the set of partitions of the set of vertices V of T . Hence, one has 

T°[T(Y^)]= r^[T%Y^)\xM^), (3.4) 

where is the Mobius function of the finite poset V{V) (see |21| ). 

2. General setting: Given a map r : C/cyc(x,x*) — > C, we define its injective version by: 
for all cyclic * -test graph T with set of vertices denoted by V, 

AT]= J2 r[T]x^,y{n), (3.5) 
7Tev{v) 

so that 

r[T]= (3.6) 

rreViV) 

This definition is then consistent with the definition of r]^ . 



torx 




Figure 5: Left: a *-test graph T (la- 
bels are replace by colors for simplicity) 
and a partition tt of its vertices, repre- 
sented by dashed lines joining vertices 
in a same block. Right: the *-test graph 



The two first assumptions we have stated in the introduction can be reformulated in terms of the 
non-injective trace. 

Assumption 1 (Convergence in distribution of traflics). 

The family of N by N random matrices Y^r has a (mean) limiting distribution of traffics r on 
Qcyci^,^*) , that is: the entries o/Yjv have finite moments of any order and for any cyclic * -test 
graph T , 



t[T] 



Hm E 



tn[T{Y 



N) 



exists. 



(3.7) 



Equivalently, by Proposition/ Definition 3.2, one can replace the trace by the injective one. 



Assumption 2 (Concentration). 

random mi 

T^[Ti(Yjv)] ...rw[T„(Y^)] 



The family of N by N random matrices Yn satisfies for any cyclic * -test graphs Ti, . . . ,T„, 



E 



N- 



t[Ti]...t[T„]. 



(3.8) 



Equivalently, by Proposition/Definition 3.2, one can replace the trace by the injective one. 



It is important to have in mind that the * -distribution of a family Yjv can be written explicitly in 
terms of its distribution of traffics T i— >■ r^r [T] . 
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Proposition 3.3 (Properties of the convergence in distribution of traffics on t/cj/c(x, x*)). 

1. If almost surely the matrices o/ Y^r are uniformly bounded in operator norm, then up to a 
subsequence Yn has a limiting distribution of traffics. 

2. The convergence in distribution of traffics on Qcyci'^,^*) of a family Yn of matrices implies 
its convergence in * -distribution. 

3. It also implies the convergence in distribution of traffics on t/cyc(x, x*) of any family of * - 
polynomials inY^. 

The first point is a consequence of a result of Mingo and Speicher recalled in Theorem |3 . 5 1 below . 
If Yjv has a limiting distribution of traffics r, the third point of the Proposition gives a sense of 
t[T], where T is a *-test graph whose edges are labelled by polynomials in the indeterminates. 

Proof of 2. Let P be a *-monomial of the form 

where iiT > 0, 7 : {1, . . . , {1, . . . ,p} and e : {1, . . . , — > {1, *}. Let Tp be the cyclic 

*-test graph whose vertices are 1, . . . , if and whose edges are (1, 2), . . . , (if — 1, if), (if, 1), the 
edge (A:, fc + 1) being labelled x^!^^^-^ (with indices modulo K). Then, one has 

^Tr[P(Y^.)] -r^[Tp(Y^)], 

so that the * -distribution of Y^r on *-monomials is the restriction of the distribution of traffics of 
Yjv on a subset of cyclic *-test graphs. □ 

Proof of 3. Let Pi,...,Pq be *-monomials and define Zn ~ (Pi(Yjv), ■ • ■ , -P,j(Yjv))- For any 
n = 1, . . . , we write 

" ^ 7,i(l) ' ■ ' 7r,(-fS'„)' 

where > 0, 7„ : {1, . . . , K^} ^ {1, . . . ,p} and e„ : {1, . • . , K^} {1, *}. 

Then, for any *-test graph T, one has TN[T{Z]y)\ — T]y[T{Y]y)\, where T is the cyclic *-test 
graph obtained from T by replacing for any n = 1, . . . ,q and e in {1, *}, the edges labelled 
by a consecutive sequence of edges ei, . . . , eK„, where Cj is labelled {^y"^('f)Y 1 j = li ■ • ■ i Kn (with 
the convention (a:*)* — Xj). The convergence of Yjv implies the convergence of Z^r. We get the 
expected result by multi-linearity. □ 



3.1.2 Statement of Assumption 3 

Recall the definitions of Mingo and Speicher [TB] . 

Definition 3.4 (Tree of two-edges connected components of a *-test graph). 

1. A cutting edge of a * -test graph is an edge whose removal would result into disconnected 
subgraphs. A two-edge connected * -test graph is a * -test graph without cutting edges. A two- 
edge connected component of a * -test graph is a subgraph which is two-edge connected and 
cannot be enlarged to a bigger two-edge connected subgraph. 

2. Let T a * -test graph. Its tree of two-edge connected components T(T) is the directed graph 
defined as follow. The vertices of T(T) consists in the two-edge connected components of 
T. Two distinct vertices of %[T) are connected by an edge if there is a cutting edge between 
vertices from the two corresponding components in T. Hence T(T) is always a tree, i.e. a 
connected graph without cycles. 

3. A tree is trivial if it consists in only one vertex. A leaf of a non-trivial tree is a vertex which 
meets only one edge. By convention, we say that the trivial tree has two leafs. 
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4-. For any * -test graph T, we denote by t(T) the number of leaves of %[T). 

Mingo and Speicher have proved in |18| an optimal estimate, reformulated in the language of *-test 
graphs as follow. 

Theorem 3.5 (Sharp bounds for the trace of test graphs in matrices, |18|1. 

Let T be a *-test graph. Let 1{T) be its tree of two-edge connected components and denote by r(T) 
its number of leaves. Then, for any family Yn of N by N matrices, 



r^[nY^)] ^iV^f^'/^-^nil^SlI' (3.9) 



where \\ ■ \\ stands for the operator norm. Moreover, there exists matrices for which this bound is 
reached. 

Their result sheds light on the quantity r(T) which turns out to plays an important role in the 
asymptotic traffic-freeness theorem we prove in this article. 

Assumption 3 (Control of growth). 

The family of N by N random matrices Yn satisfies: for any (non cyclic) * -test graphs Ti, . . . , T„, 
there exists a constant C such that 



E 



tm[T,{Ym)\ . . . rr,[T^{YM)\\ - o{N<^^y^~^ . ..N^^-^^-^). (3.10) 



Equivalently, by Proposition/Definition 3.2, one can replace the trace by the injective one. 
3.1.3 Proof of Proposition 2.2 

We go back to the proof of Proposition |2.2| which tells situations where assumptions are satisfied. 
Proof of 1 . This claim is a consequence of Mingo and Speicher Theorem |3.5[ □ 



Proof of 2. We just remark that for any *-test graph T there exists a polynomial P such that 
rjv[T(DAr)] = iTr[P(DAr)]. □ 

Proof of 3. By [TS], the convergence in distribution of traffics of Yjv and the weak local convergence 
of a Gat are equivalent. Moreover, for any *-test graph T, there exists a cyclic *-test graph T such 
that rjv[T(GAr)] = tjv [T(GAr)] . Hence, Assumption 2 impHes Assumption 3. □ 

3.2 The trafRc-asymptotic freeness of large random matrices 
3.2.1 Definition 

Recall the definition from [TB] (in an slight different formulation). 

Definition 3.6 (Traffic-asymptotic freeness). 

1. Free product of *-test graphs: Let Xi, . . . ,Xp be families of variables. A * -test graphs T 
with labels in x = (xi, . . . , Xp) is said to be a free product whenever it has the 

following structure (see Figure^. Denote by Ti, . . . ,Tk the connected components ofT that 
are labelled with variables in a same family. Consider the undirected graph Gred{T) defined 
by: 

• the vertices of Gred{T) are Ti, . . . , Tk with in addition the vertices vi, . . . ,vl of T that 
are common to many components Ti, . . . , Tk, 

• there is an edge between Ti and Vj if Vj is a vertex of Ti, i — \^ . . . , K , j — 1, . . . , L. 
Then, T is a free product m Xi, . . . , Xp whenever Gred{T) is a tree. 
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2. TrafRc-asymptotic freeness: Let Xi,...,Xp be families of N by N random matrices, 
whose entries have all their moments, and having jointly a mean limiting distribution of 
traffics on Qcyc{^,^*) , that is: for any cyclic * -test graphs T in variables Xi, . . . ,Xp, 



t[T] := lim E tat [r(Xi, . . . , Xp)] 



exists. 



We say that Xi, . . . , Xp are asymptotically traffic-free whenever: for any cyclic * -test graphs 
T in variables xi, . . . ,Xp; 

• if T is a free product in Xi, . . . , Xp, then 

f 

where the product is over the connected components ofT that are labelled with variables 
in a same family. 



otherwise, t°[T] = 0. 



Figure 6: Left: a *-test 
graph T free product in 
three famihes, represented 
by different colors. The 
non trivial connected com- 
ponents labelled by a same 
family of variables are en- 
circled. Right: the graph 
Gred{T). For convenience, 
the vertices corresponding 
to components of T are 
marked with a dot of the 
corresponding color. 



3.2.2 A traffic-asymptotic freeness theorem on tJci/c(x,x 



cyc\ 

Theorem 3.7 (The asymptotic freeness of X^^'', . . . ,Xp^-* on C/cyc(x, x*)). 

Let X^^\ . . . jXp^'* be families of N by N random matrices. Assume the following. 

1. Statistical independence: 

The families X^^\ . . . , Xp^-* are statistically independent. 

2. Joint invariance by permutation: 

For any permutation matrix Un, and any j = 1, . . . ,p except possibly one, 

[/^xf)f/i^4xW. (3.11) 

3. Convergence in mean distribution of traffics on C7cyc(x, x*): 
For any j = 1, . . . ,p, the family X^^^ satisfies AssumptionlA 

4-. Technical condition: 

For any j = 1, . . . ,p, the family X^^'' satisfies Assumptions 2 andj^j 

Then, the joint family . . . jXp^"*) has a mean limiting distribution of traffics on Qcyc{^,^*) . 

The families of matrices 'X\^\ . . . , Xp^'' are asymptotically traffic-free. 
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This theorem is a sHght different version of [16, Theorem 12.1], where we do not assume Assumption 
3 but the convergence of the distribution of traffics for all *-test graphs. 

Proof. Let T — {V, E, 7, e) be a cyclic *-test graph in variables Xi, . . . , Xp. For i ~ 1, . . . ,p, we 
denote by Ti^k, k — 1, . . . , Ki the connected components of T labelled by variables in x^ that are 
not reduced to a single vertex. In general, these *-test graphs are not cyclic. For i = 1, . . . ,p, 
denote by Vi the union of the vertices of T^ fe ioi k = 1, Ki. Then, by p — 1 applications of |16l 
Proposition 12.3], one has 



E 



(^-1)! ^ iN~\V,\)l...iN-\Vp\)\ ^^^^^_, 



iN-\v\y. {N -iy....{N -ly. 



E[n<[ri..(xl^')]] x...xE[n<[Tp,.(xW)]]. (3.12) 



fe=i 



k=l 



If T is a free product of *-test graphs in variables xi, . . . , X2, then Ti^k are cyclic for any i — I, . . . ,p 
and k — 1, . . . , Ki. By the convergence in distribution of traffics on Gcyci'^, x*) of Xi, . . . , Xp sep- 
arately (and Assumption [1]) , we get the convergence of E[T^[r(XAr)]] to the expected limit with 
minor modification of the proof of |16| Theorem 12.1]. 

From now, we assume that T is not a free product and show the following (it will be useful 
later), with the same notations as above. 

Lemma 3.8 (Tightness on the setting of the asymptotic traffic-freeness Theorem on C/cyc(x, x*)). 
With "Kpf as in Theorem \3. 7| except that Assumption^ is not satisfied, for any cyclic * -test graph 
T which is not a free product, the quantity E[r]^ [T(X7v)] | is 0[N^^) as N goes to infinity. 

First, remark that 



{N-iy {N-\v,\y...{N-\Vp\y^^^^, 



{N-\v\y 



[N -iy....{N ~iy 



. N p 



A' -1 



N 



K 



^ + ...Kp + \V\-{\V^\ + ■■■ + \Vp\)-l ^. ^-p_ 



Write \V\ = vi + W2, where Vi is the number of vertices of T that belong to a single connected 
component Ti^k, i — 1, ■ . ■ ,p and k = 1, . • ■ , Ki. Similarly, for any i = 1, . . . ,p and k ~ 1, . . . , Ki, 
write \Vi\ — v[ 



(0 



where v{ is the number of vertices of Ti^k that do not belong to other 



connected components. The number of vertices of Gred(T) is u^ed — Ki 
number of edges is Cred = k '^2'^- Hence, 

-p := Ki + ...Kp + \V\~{}Vi\ + --- + \V.p\)-l 



+ Kp + V2 , and its 



Ki + ...Kp + vi+V2-^{v[ 



- 1 



i.k 



■Kp + V2-^ V^2^ - 1 = Vred - dred " 1- 



i.k 



By the relation between the number of vertices and edges in a connected graph |12l Lemma 1.1], 
p is the number of cycles of Gred{T), that is the maximal number of edges than can be removed 
from Gred{T) while keeping a connected graph. 



On the other hand, by Lemma 3.2 for any i = 1, . . . ,p and k = 1, . . . ,Ki, one has the relation 

<[T,,,(XW)]= r^[T^^(xf^)]Mv.,(^), 

TTev(v,,k) 

where Vi^k stands for the set of vertices of Ti_k. Since for any tt in VlVi^t) one has 'c{T^^) ^ '^(Ti.fc) 
by Assumption |3] there exists a constant C such that 



E 



r^[T,,(Xf))] 
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So, by the formula (3.121 for E[T^[r(X7v)]] and the equivalent of the normalizing factor, we get 
that 

E[tO [r(XAr)]] - 0{N-P) X 0(Ari:.,.MT.,.)/2-i)^ ^. o{N-^). 

Since T is not a free product, there exists a cycle on Gred{T). Moreover, T being cyclic, each *-test 
graph Ti^k whose tree of two-edge connected components has £ leaves is responsible of the addition 
of (£ — 2)/2 cycles in Gred{T), so that the total number p of cycles in Gred{T) satisfies 



1+E 



c(r,,fe) - 2 



i.k 



Hence we get that 5 > 1, so E[t^ [T(Xjv)]] = 0(iV"i) as expected. 



□ 



3.2.3 The limiting distribution of traffics of a single heavy Wigner matrix 



In this section we prove Proposition |2.4| and show that a single random matrix satisfies Assump- 
tions 2 and 3. 



As we consider a single Hermitian matrix, it is sufficient to consider *-test graph of the form 
T ~ {V, E) (the maps 7 and e are trivial). By invariance of by conjugacy by a permutation 
matrix, 



E 



1 

N 



4>:V^{l,...,N} 
injcctivc 



J2 E Y[XNiHe)) 



eeE 



{N-\V\)l 



xS%[T{Xn)], 



where (55^[r(Xjv)] = IE[nee£; ^Jv(</'(e))] does not depend on the injective map <j). For any 1, 
denote by the number of vertices of T where are attached k loops. For any fci ^ li denote 
by (Zfei.fca the number of pairs of vertices with fci edges attaching these vertices in one way and ^2 
others in the opposite direction. Then, by independence of the entries of Xj^, 



5%[T{Xm)\^\{ 



k>l 



7V5 



n 



9fci,fc2 



Denote 



which is the number of egdes of T when the multiplicity and the orientation are forgotten. Then, 
one has 

j^OnH {Xn)\ = [[ 



k>l 



N2- 



N2- 



Since the entries of Xn are centered, 6^ [T{Mn)\ vanishes as soon as an edge of T is of multiplicity 
one, orientation forgotten. By the Cauchy-Schwarz inequality and by assumptions (1.1 1, (1.2), for 
any fc ^ 1, 



= 0(1) 
= 0(1). 
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Hence, we get that -^(5^^ [T(X7v)] is bounded. Moreover, if T has no loops and all its edges are 
of even multiplicity, then 



n- 



(3.13) 



Recall that B is the number of edges of T when multiplicity and orientation of edges are forgotten. 
By the relation between number of edges and vertices in a connected graph |12l Lemma 1.1], 

B + 1 

with equality if and only if the graph obtained from T when we forget the multiplicity and the 
orientation of its edge is a tree. In that case, we say that T is a fat tree. We deduce from the 
identity 



E 



N) 



V\^{B + 1) 



oil)) X —S%[TiX^)] 



that E[r]^ [T(X7v)]] is always bounded. If T is cyclic, since cyclic fat trees have even multiplicity 
of edges, we get by ( 3.13| 



E 



<32fc -1 

''h is a fat tree 



0(1) AT] 



It remains that satisfies Assumptions 2 and 3. Let Ti, . . . ,Tk be *-test graphs. Let T — {V, E) 
be the *-graph (*-test graph without the connectedness condition, for which trace and injective 
trace are defined by the same formulas) obtained as the disjoint union of Ti, . . . ,Tk- By |l6| 
Lemma 11.7], 

t%[T,{Xj,)\...t%%{Xm)\ = E]V^<[^"(^^)]' 



where the sum is over all partitions tt on that contain at most one vertex of each Tk, k — 1, . . . ,n. 
For any such a partition tt, denote by T^,...,T^^ the connected components of T'^ . By the 
independence of the entries of Xn, 



E 



r%[T,{XM)]...r"N[TniXN)] 



-E 



r°[r-(X^)] ...E TO[r-jx^)] 



Each expectation is bounded and converges as N goes to infinity if the *-test graphs are cyclic. 



We always has rriTr ^ n, expect for the trivial partition. Hence, E[t^ [ri(XAr)] . . 
bounded, and if the * -graphs are cyclic we get 



E 



rj^[Ti(X^,)]...<[r„(X^)] 



N- 



r°[ri]...T°[T„]. 



'N 



[T„{Xn)]] i 



IS 



3.3 Some consequences 
Corollary 3.9 (The Wigner case). 

Consider a family X^r of independent Wigner matrices, independent o/ Yjy. Assume that 
converges in * -distribution and that it satisfies Assumptions 2 and 3. Then, (Xjv, Y^v) converges 
in * -distribution and are asymptotically * -free. 

Proof. By Assumptions 2 and 3, Y^r is tight on tJcj/c(x, x*). Consider a subsequence along which 
Yjv converges. Let X^ be a heavy Wigner matrix with trivial parameter (a, 0, 0, ... ). By Theorem 



2.1| (Xtv, Yjv) has a limiting * -distribution along this subsequence, given by formula (4.5 1. The 
family X^v converges to a family of semicircular traffics, traffic free from the limit of Y^y. Since 
the traffic-freeness of semicircular variables with arbitrary traffics implies their *-freeness [17[, we 
get that Xat and Yjv are asymptotically *-free variables. □ 
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Proposition 3.10 (Heavy covariance matrices). 

Let N be an integer. Let Nq^\...,N^^ he integers such that N^j^^ ^ CkN, Ck > for any 
k — 0, . . . , K . Let Wtv = {Wi'^'^ , ■ ■ ■ , W^^"*) be a family of random matrices such that: 

• for any j ^ 1, . . . ,p, one has VF)^^ = Mj^'zj^^il/j^^* . 

• Mjv is a family of independent random matrices AI^^\ j — 1, . . . ,p, with independent entries 
having the same distribution of a random variable rrij^N such that 

for any n ^ I. The matrix Afj^' is of size Nq by Nj^^^ for an integer kj in {1, . . . ,K}, 
j = 1, . . . ,p. 

• Zj"'^'' — (z[^\ . . . , Zp'^^) is a family of random matrices, and the matrix z'^^^ is of size Njf^"^ 



by Tvf ) 



Let~YN be a family of Nq by Nq random matrices. Assume that 

1. the families of matrices Mpj, Yjy, (-^j"'^^)fej=fe- k = 1, . . . , K are independent, 

2. the families of matrices Yjv, {Z^'^'^)kj=k, k = 1, . . . , K , satisfies the assumption of Theorem 



\2.1\ separately. 

Then, the family of matrices (Wjv, Yjv) has a limiting * -distribution as N goes to infinity. 

Proof. We prove the Proposition for K ~ \, the result can be obtained by recurrence on the 
number of size of matrices. Consider the square matrices of size [Nq + iVi), by blocks 

We consider the matrices 



eo= ^° ^ ], ei= 1.. h ^, - 'an)* ^ln) ,J = 1,-.-,P, 



where Xj^'^\Xj^'^\ j = are independent heavy Wigner matrix with parameter 

( lim E[7VTO2V])fc^i> 

square of size Nq and N\ respectively, independent of (Yat, Zat, Mjv). It 

can be noted that 

■rjV [T(Yjv, Zat, ei, 62)] = TaTq [r(Yjv, 62)] X -^——Ij^ labelled in (y.es) 

ci ' 

+ tatq [T(Zjv, ei)] X c It labelled in (z,ei)- 

Co 



Hence, (Y^r, Z^r, ei, 62) satisfies the assumption of Theorem 2.1 and so (Xat, Yat, Zjv, ei, 62) has 
a limiting * -distribution, given by the traffic free product. Since for j — 1, . . . ,p. one has Mj'^^ = 

eoXj^^ei and W"j^^ = Mj^^zj^^Afj^^*, we get that (WAr,YAr) has a Hmiting *-distribution. 
Moreover, for any *-polynomial P 

TrfF(Wjv, Yjv)] = ^ ^TrfPfWAr, Yat)] . 

Ni+No ^ ^ J^' J^U i^ + l^o ^ ^ 

Hence the convergence of (Wat, Yat). □ 
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Proposition 3.11 (Independent Levy and random matrices). 

Let Xjv = {x[^\ . . . , Xp^') be a family of independent Levy matrices, independent of a family of 
random matrices Yjv satisfying the assumption of Theorem \2.1\ Then, for any Hermitian matrix 
Hf4 = P(XAr,YAr), where P is a fixed * -polynomial, the empirical eigenvalues distribution of Hn 

converges weakly, i.e. for any continue bounded function / : R ^ R, E -^Tr[/(iJAr)] converges. 



Proof. Let _B ^ 1 be a positive and large number, and set for any j — 1, 



X 



(B.N) _ 



(X]^^(m,n) X l|^(iv) 



{m.n)\^B J m,n— 1, 



,N' 



1.2.1 



the family 



is a family of independent heavy Wigner 



Y^f) has a limiting * -distribution. Hence, the empirical 



By Section 

matrices. Hence, by Theorem 2.1 

{ E T — 

eigenvalues distribution of Hj^ = 

( B ) 

reasoning as in [5 , Section 8], Hj^ ' and i/^v are closed is the sense of rank and by [5, Lemma 12.2], 
this implies that the empirical eigenvalues distribution of Hj^ converges weakly to a measure /i 
and /J, = limB_j.oo □ 



P(X^\ Yat) converges weakly to a measure /i^^-*. By the same 



4 Limiting ^-moments of independent heavy Wigner and ran- 
dom matrices 

4.1 Proof of Proposition 2.8 



Let Xjv = {X[^\. . . , X^^'>) and Yn be as in Theorem 2.1 We denote by $ their mean limiting 
* -distribution and by r their limiting distribution of traffics on 5cj/c(x,x*): for any *-polynomial 
P and any *-test graph, 

$(P) = lim lTr[P(XAr,Y^)], (4.1) 

t[P] = lim TAr[r(XAr,YAr)]. (4.2) 

The family of variables x — (xi, . . . ,Xp) corresponds to X^r, the family y correspond to Y^r. 
Theorem 2.1 tells us how $ can be written in term of the injective version of r, defined by (3.51. 
Consider a polynomial P of the form 

P^x^^^P,{y)...x^^L)PL{y), (4.3) 

where 7 : {1, . . . , i} — > {1, . . . ,p}. Let Tp be the *-test graph in variables x, Pi(y), . . . , Pl{y) as 
in Figure |7] 



Figure 7: Example with L = 4. Starting from 
~ the vertex marked with a dot and running clock- 

wisely, the edges are labelled a;^(i), Pi(y), . . . , x.^^^ 
and Phiy)- Then, the black edges have la- 



( ) 

bels a;^(i), . . . , and the red ones have label 

C> Pi(y),...,Pi(y). 

Then, by the trafhc-asymptotic freeness of xj^'', . . . .x'^^ and Y^r, we obtain the formula 

^{Tp is a free product) nr"[f], (4.4) 
ne'P{2L) f 



where 

• P(2L) denotes the set of partitions of {1, . . . , 2L}, 
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• Tp is the *-test graph defined from Tp and tt as in Figure [5] 

• the notion of free product of *-test graphs, given in Definition |3.6[ is relatively to the family 
of variable (a;i), . . . , {xp), (Pi(y), . . . ,Fft:(y)), 



t' is the injective version of r, defined by (3.5 1 



• the product is over all connected components of Tp that are labelled by a family among 
(xi),..., {xp), (Pi(y), . . . , ^L(y)), as illustrated in Figure [ej 

We have drawn two examples of free products of *-test graphs in Figure|8] remembering the marked 
point of Tp as we did in Figure [7] 

Let TT be partition in V{2L) such that Tp is a free product. Let be the *-test graph ob- 
tained from Tp by merging the components labelled Pi(y), . . . , Ptly)- Hence, must be a fat 
tree. From tt and F'^ , we get a minimal cycle on a tree (Gtt, c^) in as we take care of the way 
we fold Tp into Tp. The two partitions of Figure [8] give the same minimal cycle. 





Figure 8: Two free products of fat trees (black) and arbitrary *-test graphs (red) for a polynomial 
P ^ 2:^(1) Pi (y) . ..x^(L)PL{y) with L = 8. 



Let (G, c) in L'''^\ All partitions tt such that {G-,^,Ct;) = {G,c) will give the same contribution 
ujHwiG, c) from heavy Wigner matrices that can be factorized in Y[f ''''^i'^]- 

$(P)- J2 ™(G,C) ^ l(G,,c„) = (G,c)n^°[^']' 
(G,c)e£<^) 7veV{2L) T' 

where the product on T' is now over all connected components of Tp that are labelled Pi (y ),..., (y) . 
Let (G,c) in L'^'^K It remains to show that 

l(G.,,,) = (G,c)n^"[^'] =^™(G'C). 

■KeV(2L) T' 

For any v vertex of G, recall that we have defined a *-test graphs Ty labelled in Pi(y), . . . , Pff(y). 
All the partitions tt such that (Gt^Ctt) — {G,c) give the same fat tree F — F^ . "Replace" the 
vertices of F by corresponding *-test graphs T„'s in the following way: 

1. consider the disjoint union of the Ti,'s. 

2. By construction, each vertex of a r„ is associated to an edge of G. Link the vertices of two 
different T^ and that correspond to a same edge of G by n edges, where n is the number 
of times c walks on this edge. 

3. Orient half of these edges in one direction and the other and the other direction. 
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4. Label these edges by the color of the corresponding step of c. 

The *-test graph we obtain is Tp° , where ttq is the coarser partition for which {Gjr„,CT^„) = (G,c) 
(see the rightmost *-test graph in Figure [s]). The other partitions tt which give (GTrjCTr) = (G,c) 
are the sub-partitions of ttq which do not put in a same block vertices from different T„'s (as for 
the leftmost *-test graph in Figure |8] compared to the rightmost). By the relation ( |3.6| between 
the standard and the injective trace, we obtain the expected result and hence the Proposition. 

4.2 Proof of Corollary 2.9 

Proof of 1. and 2. As a cycle visits a tree with different colors for each vertices, the total number 
of steps of a given color is an even number. We then obtain that ^{xj) = for any j = 1, . . . ,p 
and the second point of the corollary. Moreover, <i>(a;|) is Oj.i since there is only one cycle running 
on a tree with one edge in two steps, which gives the contribution aj^i. 

Proof of 3. The set of cycles running on a tree with rii steps of colors zi, then n2 steps of 
colors 12, and so on, is in bijection with the product of the sets of cycles running on a tree with rij 
steps of colors ij, j — 1, . . . , L as soon as the colors are distinct. The weights ujhw for x^^ . . . x^^ 
are the products of weights for x"^ ,. . . ,x"^ . The weights ujtr comes from the root of the trees. 

□ 

4.3 The false freeness property 

The folding trick gives a unformal algorithm for the enumeration of the cycles running on trees 
when L is not to large, and then for the computation of limiting joint moments in heavy Wigner 
and deterministic matrices of few degree. Denote by ^ the elements of such that the cycles 
visits each edge of their tree twice. They correspond to cycles that contribute for classical Wigner 
matrices, a case where Xjv and Yjv are asymptotically free in the sense of Voiculescu [21_. The 
elements of \ £2^' are enumerated by folding the combinatorial objects of £2^''. 




Figure 9: Left: a couple (G, c) and its unfolded version. Right: starting with a double graph, all 
the couples (G, c) obtained by the folding trick. 



Unfolding trick: Let (G, c) in £W After some steps, leaving a vertex v, the cycle c comes 

back in an edge it has already visited. Then it induces a sub-cycle c on the tree of the descendent 
of V. We create a copy G of the sub-tree induces by c, forget its original embedding and embed it 
in such a way c respects the rules concerning the order of visits of the edges of G. Then we attach 
G endowed with this new orientation at the vertex v, between the edges it has already visited and 
the others. If some edges of the tree of the descendent of s where only visited by c, then we erase 
them. We then keep an element of C^^\ Iterating this procedure a finite number of times, we then 
get an element of £2 ^ . 

Folding trick: Reciprocally, let (G, c) be an element of C^^K Chose an edge ei of the tree. 
If possible, chose an other edge 62, which shares the same vertex toward the root and which is of 
the same color as ei. Then, merge these two edges, draw the tree of the descendant of ei at the 
right of the the tree of the descendant of 62 and redirect the cycle c in this new tree. We then 
obtain an new element of £^^^ For any element (Go, cq) of £3^^ we denote by fold{Go, cq) the set 
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of all elements of we get by applying many times this trick. 

Folding and unfolding tricks are illustrated in Figure [o] Two different elements of £2^-' have different 
folding sets. We then get from this construction the following proposition. 

Proposition 4.1 (The false freeness property). 

For any * -polynomial P of the form P = a:^7(i)-Pi(y) . . • a;7(L)^'L(y), one has 



E E ujhw{G,c)xu:tr{G,c), 

(Go,co)e£<"'' (G,c)e/oid(Go,co) 



(4.5) 



where ojhw md ujtr o,re as in Proposition^ 

The false freeness property gives a method to reasonably compute limiting joint moments of heavy 
Wigner and deterministic matrices: 

1. Enumerate the elements of €'"2'^. 

2. Fold the branches of these colored trees. 

3. Then, read the contribution of each element. 

We do not describe how to be sure to obtain all the elements of C'^'^^ during the second step of the 
algorithm, as our purpose is to use this method for relatively small L. As an example, we have 
computed ) = 5af + 601^201 1 + 01,3 in Figure 



10 





Figure 10: Computation of : ^{xf) — 5af i + 601^211 1 + ai,3- In the frame, we have enumerate 
the elements of c'^'"''^\ Three of them, the rightmost ones, can be folded in an unique way. The 
double tree on the bottom of the frame can we folded in four ways, as in Figure [9j On the left 
of each minimal cycle on a tree (G, c), we have drawn the permutation tt of {1, . . . , 6} such that 
{GTr,CTr) = (G, c), with the notations of the proof of Proposition 2.8 It should be noticed that 
the partitions on the frame are the dual, in the sense of planar partition, of non crossing pair 
partitions. Planarity is broken when trees are folded. Then, these partitions are the dual clickable 
partitions (see |22|). 
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4.4 Proof of Propositions 2.10 and 2.11: the non asymptotic *-freeness 
of Xjv and Y^v 



Proof of f{xi, X2) = ai,2a2.2; where xi,X2 are heavy Wigner of parameters (ai,fc)fe^i, (a2.fe)A;^i . We 
first expand the quantity 

- <^{x\)) {xl - ^xj)) {xl - HA)) {4 - H4))) 

= <^{xlxlx\xl) 

^2Mxl)'^{xlxi) - 2<^{xl)'^{xlx\) 
+4$(:r2)^(^2)^(^2^2) ^ $(2:2)2$(^4) + Hx2fHx\) 

-A^xlf^xlf 
+Hxl)H{xlf. 

Using the traciality of $ (that is ^{PQ) — ^{QP)) and interchanging the roles played by xi 
and X2, it is enough to compute ^{x\), ^{x\x\). $(a;fa;|) and ^{x\x\x\x^). By the basic 

we have <^{x\) — ai_i, ^{x-fx-^) = $(a;i)<i>(x|) — 01^102,1 and 



2.9 



properties stated in Corollary 




Figure 11: Computation 
of ^(xf) = 2a?_i +ai,2. 



Figure 12: Computation of ^(x^Xgxf = Saf lO^ ^ + 01,202 1 + "^i i'^2,2 + ai,2a2,2- Edges corre- 
sponding to Xl are in black, the ones corresponding to X2 are in red. For technical reasons, we 
have drawn fat trees instead of cycles running on trees. 



The computation of ^{xj) and ^{xfx2XiX2) are done in Figures 11 and 12 following the algorithm 



of the false freeness property (Proposition |4.1[ ). As we do not consider deterministic matrices, we 
do not need to open boxes in the vertices as in Figure [3j This gives 

- $(.?)) (xj - HxD) {xl - HxD) {xl - HxD) 

~ 3a^ ;^a2_;^ + 01,202.1 + Oj 102,2 + 01,202,2 
-2ai,i(2a2 1 + 02,2) - 2a2,i(2ai ^ + 01,2) 
+4ai,i02,i + Oi,i(2a2 1 + 02,2) + a2,i(2ai ^ + 01^2) 

q 2 2 

'^1,1*^2,1 
= 01,202, 2- 
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□ 



Proof of Proposition \2.1C\ and Proposition \2.11\ 1 . We expand the quantity 
ci>((x? - (y^ - $(y^)) {xl - {y' - ^(y^))) 

We have to compute ^(xi ), $(a;i), ^{xfy"^) and '^{x'^y'^x^y'^). Using again the basic 



properties of Corollary 2.9 and the computation of *&(a;^) of Figure [ll] the only term we have to 
compute is ^{xfy^xly"^). 




g 



8 





Figure 13: Computation of $(xiPi(y) . . . xiP4(y)) . We consider the enumeration of cycles of 
length 4 running on a tree of Figure [TTj and them open boxes on vertices. 



In Figure 13 we have computed $(a;iPi(y) . . . a;iP4(y)) for any *-polynomials Pi, 
obtained 

<i>(a;iPi(y)...a;iP4(y)) = a?,i<i>(Pi(y)P3(y))<I>(P2(y))$(P4(y)) 

+a?,i$ (Pi (y)) $ (P2 (y)P4 (y)) $ (Pa (y)) 
+ai.2$('n^i(y)'^3(y))$^'n^2(y),P4(y)). 

Specifying Pi = P3 1 and P2 — P4 ^ y^, we get 



^xjy'xly') = a\ Mv'f + a? Mv') + a^^^'Hy', 2/') 



We then obtain 



((.? - (y^ - ci>(y2)) {x\ - HxD) (y^ - My')) 

= a?,i$(y2f +a?,,$(y4)+al,2<i>(^)(2;^y2) 
-2a?,i$(y4)^2<i>(2;2)2(2a2^+ai,2) 
+4a2 i$(y2)2 + i$(y4) + $(y2)2(2a^_^ + a,^^) 

-3alMyr 

= ai,2($(2)(y',y') - Hy'?) - ai.2<i>(^n2/' - - 

where in the last equality we have used the bi-linearity of 



, P4 and 



□ 



Proof of Proposition \2.10\ 2. If Y/v is a diagonal matrix, then it turns out that <i>(^) (P, Q) = <^{PQ) 
for any *-polynomials P, Q, and the claim follows directly. □ 
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5 Proof of the Schwinger Dyson system of equations for heavy 
Wigner and diagonal matrices, Theorem 2.12 

The idea of the proof of Theorem |2.12| is to classify, in the combinatorial approach by cycles run- 
ning on trees, those for which the cycle visits a fixed number of times the first edge of the tree. 



Before providing Theorem 2.12 we apply theses equations to give an other computation of $[ 



— {xi X 1 X Xi)x2XiX 



First, we enumerate the decompositions 

{Xl X Xix\ X XljXlX^ 

{xi X xix\xi X Xl)X 
{xi X 1 X xi)x\{xi X 1 X xi)x 



Then, by Theorem |2. 12| we get 

^[xfxlxlxl] = ai^i(^$[l]<i>[a;2Xia;2] + <I'[a;ia;2]$[a;ia;2] + <i>[a;ia;2a;i]$[a;2] j 

+ai,2$(^)(l,l)$(^)(^i^^) 
= aia($[a:?]$[a;4]+0 + <i>[x?]$[a;2]2) +au2^^^H4,4) 

= aiiali+al^^x^]+ai,2'i>'^^\xlxj), 

where we have used the facts that $[a::"a;™] = $[a::"]<i>[a;™] for any n, m ^ 1 and that ^[xf 
for i — 1,2. By Theorem |2.12[ one has with a similar computation 

^x4] = a2s[^mxl] + <i>[x2Mx2] + 'P[x2]<i>[l]^ +a2,2$^'ni,l)**'^(l,l) 
= 2a2i+a2,2- 



To compute $(^)(a;2, -"z^i) with Theorem 2.12 we enumerate the decompositions 

(a;2,X2) = (^{x2 X 1 X 2:2)1, ^2^ = (^{x2 X 1 X .i;2)l, lix2 X 1 X 2:2)1^ 

So we have 

<f^^H4,4) = a2,i$[l]$(')(l,a;2)+a2,2$(')(l,l)$('^(l,l,l) 
= 02,1+02,2- 

We then get as expected 

^[x^xlxlxl] = ;^a2 ;^ + a^ ;^(2a2 ;^ + 02.2) + ai,2(a2,i + 02,2) 
— Sal ial i + of ia2^2 + ai,2a2,i + ai,2a2,2- 



Proof of Theorem \2.1/^ For clarity of the exposition, we start by proving (2.6 1 for ivT = 1, that is: 
for any j = 1, . . . ,p and any monomial P, 

<i>(a:,P) - ^ flj- fc (5-1) 
fc^l XjP—{xjLiXj)R-i...{xjLkXj)Rk 



We write P = Pi(y) x a;..y(2)P2(y) . . • x^i^L)PL{y) and set 7(1) ^ j. By Proposition 

<i>{XjP)^ UJhw{G,c) X UJTBiG.c), 

(G,c)e£<"') 

where ujhw a-nd ujtr are given in Definitions |2.6| and |2.7[ 

Step 1: Cycle visiting 2K times the first edge 

Let (G,c) in C'''^\ The root of G is called the vertex number 1, the second vertex visited by 
c is called the number 2. Saying that the undirected edge {1,2} is visited exactly 2K times is 
equivalent to say that 
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1. there exist cycles d'^^\ . . . , starting at the vertex 2, 

2. there exist cycles e^^-', . . . ,e*^^^ starting at the vertex 1, 

3. theses cycles do not visit {1,2}, 

4. c can be written 

c = a o o a* o e^i) o a o dP^ o a* o e^^) o • • • o a o df^^^ o a* o e^^\ (5.2) 



where o denotes the composition of paths, a = (1, 2) and a* = (2, 1). See Figure 14 for an example. 





Figure 14: Left: the cycle visits 4 times the first vertex. Right: the two couples of cycles induces. 



Assume that c is of this form. Since the edge {1, 2} can only be visited by steps of color 7(1), we 
can write 

(7(1),...,7(L)) = (7(1), aW, 7(1), 7(1), ^(1)^^(2)^. (5.3) 

where for any k = 1,...,K, a^^) = (a^^^ (1): • ■ • , and ^'^''^ = (/3('=)(1), . . . ,/3W(i^,fc)) 

are families of integers in {1, . . . ,p} for any k — 1, . . . ,K. The integer La,k is the length of the 
cycle a^*^^ and Lp^k is the length of the cycle f3'-''\ 



The two cycles d = d^^^ 



' d(^) and 



=(1) 



o e 



are induced by c on disjoint subgraphs 



Gd and Gg of G respectively. For any k ~ 1, . . . ,K, the cycle d''^-' has steps of color a^'^^ so that 
the steps of d are colored by a = {a^^\ ■ ■ ■ , a^-^^). The same holds for e with /3 = (/3(i), . . . 

Hence, by rooting Gd on 2 and Gg on 1, we get that {Gd,d) belongs to and (Ge,e) be- 

longs to C^^\ They are not typical elements of these sets, in the sense that for any k ^ 1, the 
cycle d always comes back to the root of its tree after La^i + ■ ■ ■ + La,k steps. Hence the following 
definition. 

Definition 5.1 (Chain of cycles running on trees). 

Let K ^ 1 be an integer, L = (Li, . . . , Li^) be a family of integers and 7 : {1, . . . , L} — > {1, . . . ,p}, 
where L = Li + ■ ■ ■ + Lx- We set C^"^ the set of all couples {L,c) in C^'^\ such that c can be 
written c = ci o • • • o ck, where for any k = 1, . . . , K the cycle Ck is of length Lk- 

If we denote Lq, = [La.i, ■ • ■ , La^x) and — {Lp^i, . . . , L^^x), then we get that actually {Gd, d) 



belongs to /^l"-* and that (Ge,e) belongs to C'^ 



(/3) 
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Step 2: Reciprocal construction 

Let if ^ 1 be an integer and consider a decomposition 

(7(1),...,7(L)) = (7(l),aW,7(l),/3W,7(l),a(2),^(l),/3(2),...,^(l),aW,^(l),/^^^^ (5.4) 

where for any k — 1, . . . ,K one has a^^^ is in {1, . . . and Z?*-*^^ in {1, . . . ,p}^^-'' for sequences 

of integers = {L^at ■ ■ t L^^k) E^nd = [L^i, . . . ,Lpx)- Define a = (a^^^ . . . , a^^-*) and 
/3 = (/?(!),. . . , /3(^)). Let (Gd, d) in and (Ge, e) in /:[f^\ We write d = d^i) o • • • o d(^) and 
e = e^^-* o • • • o e*-^'-* with the notation for chains of cycles (Definition 5.1 1. As we embed the graphs 
Gd and Ge, Unk their roots by an new edge, extend the cycles into a cycle Cd.e as in (5.2) and set 
the root of Gd,e to be the root of Gg, we get an element (Gd,e,Cd.e) in C'^^'^ 

As (Gd,d) and (Ge,e) are the chains of cycles running on a tree of the previous step when we 
start with {Gd,e,Cd,e)i we have proved that 

^(XjF) = ^ ^ ^ LOHw{Gd,e,Cd,e) ^ UJTR.{Gd,e,Cd,e)- (5.5) 

as in j5.3| l ^ ' 

S'iep 5.' Computation of the weights 

We have obviously uiHw{Gd,eTCd,e) — aj^K ^Hw{Gdj d) x w//vy(Ge,e). For the weight ujtr, it is 
important to take care about the dependence on the polynomials involved. We then write ujtr = 
^TR.{P) in the notation of Definition 2.7 In the case considered in that proof, the polynomial has 
been denoted XjP. Recall that 

U^TR{x,P){Gd.e,Cd^e)^ \{ t[T,], (5.6) 

V vertex of Gd^e 

where the *-graph tests are obtained by opening boxes on the vertices of Gd,e as in Figure [3[ Since 
the matrices Yjv are diagonal, we actually have a much simpler expression for this weight. First, 
remark that the diagonality implies that for any vertex v of Gd,e, 't^v] = '''\T^] where tt is the 
partition of the vertices of Ty with only one block. Hence, this quantity is equal to $( Hn -F'n(y)) ; 
where the product is over all integers n = 1, . . . ,L such that the n-th step of Cd,e is incident in v. 
For V different that the roots of Gd and Ge, this quantity is the same as if we replace Cd^e by Cd 
or Ce (depending if the vertex comes from Gd or Ge). For v the root of Gd or Ge, we have to take 
into account the steps of Cd.e on the edge {1, 2}. 



Given (c?,e) as in the sum (5.5l, we get a decomposition of the polynomial XjP: 

XjP — {xjLiXj)Ri . . . {xjLKXj)RK, 



where the position of the x/s corresponds to the position of the 7(l)'s in (5.4|, where we have 
decomposed (7(1), . . . ,7(if)). As XjP — x-y(i)Pi(y) x a;^(2)P2(y) x • • • x x^[L)PL[y), we can write 
Lk = Pi^{y)LkPj^-i{y) and Rk = Pj^{y)RkPi^^,-i{y) for any k ^ 1,...,K. The integers ife's 
corresponding to the steps where the cycle runs on (1, 2), the integers jk to the steps where it runs 
on (2,1). 

We set_Lfc = LkPj^-i{y)Pi^{y), Rk = RkPi„^i-i{y)Pj^{y) and then L Li...Lk and R = 
Ri...Rk- As the polynomials Pj^_i(y) and P^Jy) (respectively Pi^^j_i(y) and P,^(y)) con- 
tribute in the root of Gd (respectively Ge), we get 

UJTRiXjP){Gd,e,Cd,e) = WTfl [i] (Gd , Crf) X ^Tfl [P] (Ge , Ce) . (5.7) 

Step 4: Conclusion of the combinatorial decomposition 
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We have obtained that 



^{xjP) = ^ aj^K Y ( Y ^^HwiGd,Cd) X ujTRiL){Gd,Cd) 



as m \5.3\ " 



^ W//VK(Ge,Ce) X a;Tfl(i?)(G'e,Ce) |. (5.8) 



The classification in the sum over a, /3 is in correspondence with the number of way we can 
decompose XjP into 

XjP = {xjLiXj)Ri . . . {xjLKXj)RK- 
It remains to interpret the combinatorial terms in terms of the multilinear forms {'^^^^)k'^i- 



Proposition 5.2 (Proposition 2.8 continued). 

For any polynomial Pi,...,Pk of the form Pk = a;^j^(i)Pfe,i(y) . . . x^^(Lk)PLt{y), k ^ 1, one has 
$(^)(Pi,...,Fk) - u:hw{G,c)xu:tb{P){G,c), (5.9) 

(G,c)e£<7' 

where L = (Li, . . . , L^), 7 = (71, . . . , 7^) seen as a map {1, . . . , Li + • • • + i/f } — )■ {1, . . . ,p}, and 
P = Pi . . . Pk- The weights ujhw O'f^d ujtr are the same as in proposition\2.8\ 



Proof of Proposition \5.^ Based on formula (4.4 1 given by the traffic asymptotic freeness of X 
and Yat, this proposition is obtained by a minor modification of the proof of Proposition [2 



AT 



Denote by Tp^ , ■ ■ ■ , Tp^ and Tp the *-test graphs obtained as in Figure[7]for polynomials Pi , . . . , Pk 
and P respectively. Consider the *-test graph T obtained by merging the roots of Tp^ , . . . , Tp^ . 
Then, <i>(-'^)(Pl, . . . ,Pk) = t[T], where r is the limiting distribution of traffics of (X^r, Y^r, YJ^). 

Let TTg be the partition of {1, ... , 2(Li + • • • + Lk)} that puts in a same block {1, 2Li + 1, 2(Li + 
L2) + 1, . . . , 2(Li + • • • + Pa'-i) + 1} and let alone the other integers. Remark that T = Tp'. By 
(4.4 1, one has 

$(^)(Pi,...,Pa) = Y 1(T- is a free product) n^°[^] 

TTeV(2{Li + ---+LK)) f 

= 1(T]; is a free product) [^'l i 

where tt ^ ttq means that tt is a sub-partition of ttq. We find the result with the same reasoning 
as in Proposition |2.8[ as we realize that the condition tt ^ ttq exactly means that the cycle comes 
back at the root after 2Li, 2L2, . . . and 2Lk-i steps. □ 

By ( |5.8[ ) and Proposition, we get 

<i>{xjP) = Y J2 <^^''HLi,...,Lk)x<P<^'^\Ri,...,Rk). 

As the matrices Y^r are diagonal, for any polynomials Pi, ... , Pk, any polynomial Q(y) and any 
k — 1, . . . ,K , one has 

$W (Pi, ... , Pfe_i, PfeQ(y), Pfc+i, . . . , Pa) = (Pi, . . . , Pfe-i, Q(y)Pfe, Pfe+i, . . . , Pk). 
Hence we get as expected 

<^{x,P) = E ^^''\Li,...,Lk)^^''''KRu...,Rk). 

K^l X j P —{x j L ix j) B,\ . . .{x j L^Xj)R,ji 
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Step 5: the general case 

Let Pi , . . . , Pk be monomials of the form 



Pi -Pi J (y) X 2^71(2)^,2(7) •• •a;7i(ii)^i-ii(y) 

Pk = x^,(i)Pfc,i(y)...a;^,(L,)Pfe,Ljy), A: = 2,...,ii'. 



By setting 7 and L as in Proposition |5.2| we have 



$W(x,/'i,...,Pk) = I] c^w(G,c)xc^™(P)(G',c). (5.10) 

(Gx)e£(7' 

Let (G, c) in /^l'''^ ^'^i^ write c as a composition of cycles of length Li, . . . , Lk, namely c = cio. . .ck- 
Saying that c visits {1,2} exactly 2k times is equivalent to say there exist non negative integers 
Si, . . . , Sk such that 

• Sl > 1, S2, . . . ,SK ^ 0, 

• si + ■ ■ ■ + Sk — k, 

• for any m — 1, . . . , K, the cycle Cm visits a exactly 2sm times. 

Assume that for any m — 1, . . . ,K, the cycle c„i visits {1, 2} exactly 2sm times. Then we get a 
decomposition 

ci = a o d(i'i) o a* o e^^'i) o a o „ ^* ^ ^(1,2) ^...^ao S^^'^'^ o a* o e^^'''\ 

and for any m — 2, . . . , K , 

c„ = e("'°) o a o o a* o e^™'^^ oao S"''^'> o a* o e'^^^^ o • • • o a o o a* o e^^'^i). 

The only difference is that the cycles ci, . . . , c™ are not constrained to visit {1,2} during their first 



step. The rest of the proof can be written as we made for the proof of (5.11, without any new 



niceties. We the same reasoning as before, we obtain the expected result, i.e. Theorem |2.12[ □ 

6 Proof of Proposition 2.13: a characterization of the spec- 
trum of a single heavy Wigner matrix 



We manipulate truncated sums. We write computations based on Theorem |2 . 1 2 1 for the truncation 

1 ^ 1 
F^{K):=jj^ -ci>(^)(2^"^+\^"%---,2;"-) 

n—O niH \-nK—Ti 

of the formal power series associated to 

3>(^)(a;(A - x)-\ (A - x)-\ . . . , (A - x)-^). 
Remark that it is equal in the sense of formal sums to 

and then it enough to prove 
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in the convergence of power series in j. We fix an integer N ^ 1. By Theorem 2.12 we have 

1 ^ 1 

71— Q niH l-riK—n 

1 ^ 1 

Xi?E E E 



n=0 niH hnK=n l<k< "+'- 



X 



siH l-SK:=fe (r,t) 



The last sum is over all families of non negative integers 

r _ ,'^(1) „(i) JK) ^(K) ^ 

t - +(1) +(2) ,(2) AK) (K)-. 

such that 

rJ^^ + --- + rW+4i) + ...+t(i) = „i + l-2si, 

^ + • • • + r« +4'^ + ■■■+ = - 2s„ z = 2, . . . , X. 
We have used (and we will use) the notation 

The restrictions on the third and fourth sums follow from consideration on the degree on the mono- 
mials we compute. Now we interchange the order of summation of (ni, . . . , hk) and (si, . . . , sk)- 

n— ni-i l-riK—n 



E E E ^ 

l^k^H+l siA \-SK=k 2k-lt^7it^N 



X 



E EE<5'^'H-'")E*^'^''"^(-*)- 



ni-\ l-TiK—n 1 r t 

ni^2si-l 

By the sum over 1, we mean the sum over all families of non negative integers 1 = (li, ... ,1k) such 
that 

0< h s^ni + l-2si, 

< h ^ ^2 - 2S2, 



^ Ik ^ n-K — 2sjf . 
By the sum over r, we mean the sum over all families of non negative integers 

V. _ ^(1) J^) r^K). 

^ — V 1 )•■•)' si ;•■■>' 1 T ■ ■ T ' Sk n 

such that 

„(1) 



r 



if- 
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At last, by the sum over t, we mean the sum over all families of non negative integers 

t _ +(1) +(2) .(2) AK) (KU 

such that 

= ni + l-2si-/i, 

t^o^+---+€J = n2-2s2-/2, 

to^^+---+tiV = nk-2sK-lK. 
Given k,si, . . . ,S2 as in the previous formula, we set the change of variable for n, ni , . . . , uk 

m = n + 1 — 2k, 
mi = ni + 1 — 2si, 
m2 = n2 — 2s2, 

triK = riK- 2sk- 



Remark first that 

111 1 1 



Hence we get 



4E E 



1 + 1 ^"2 

\K Xn 

n=0 niH \-nK=n 



JV+l-2fe ^ 

= E E E 

l^fe^^^^ siH \-SK=k m=0 

X E E Ei^^'^(--)EA^^^'=""-^n:^*). 

"ilH |-mjf=m ii=0...)ni r t 

The sum over r is the same as before, and now last sum is over all families of non negative integers 

f _ (fW f(l) y.(2) A2) AK) JKU 

I, — yo-^ , . . . , Og^ , Lq , . . . , (,^2 , . . . , (,q , ■ • • , I^sk > ' 

such that 

4^) + . ..+4^) = rUK-lK. 

We replace the set variables (mi, . . . , rriKji, ■ ■ ■ ,1k) by variables pi, . . . ,Pk and qi,. . . ,qK where 
for any i = 1, . . . , K we have set pi = rrii — U and qi — k. Then we get 

1 ^ 1 



XK ^ ^ 

n— niH \-nK=n 



N+l-2k ^ 

E E E E 

l^fe^iv±i m=0 (p,q) si + ...+s^=ft 

si^l, S2,.-->«-R:^0 

xEi^^'=H-^)EA^^^'=""-^H-*). 

r t 
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The sum over (p,q) is the sum over all families of non negative integers p = {pi,. . . ,pk) and 
q = (gi, . . . , qx) such that 

Pi H \-pK + qi + ■■■qK = m. 

The sum over r is the sum over all families of non negative integers 

r - (-^(1) ^(1) JK)^ 
'■ — V 1 I ■ ■ • I ' Si ; • ■ ■ ) ' 1 ; • ■ • I ' sjf ;i 



such that 



The sum over t is the sum over all families of non negative integers 

~ V-l I ■ ■ ■ I ''Si : ''0 1 ■ • ■ : '-S2 ' ■ ■ ■ ' 1 ■ • ■ J ''sjf /Ij 



such that 

Let iiT ^ 1 and fc ^ 1 be integers. Then there exist tuples of non negative integers 

(si, . . . , sk) such that si + . . . sk — k, si ^ 1 and S2, . . . , s/f ^ 0. Hence we get 

n — n. 1 H 1- n ^ ^ n 

fK + k-2 

2^ afe 



K -I 



0^p+9^JV+l-2fc r-iH hrfc=g 

tlH htfc + K-l=P 

This gives the expected result by identification of the coefficients. The uniqueness of the solution 
of the equations follows directly from the observation of the valence of the formal power series. 

A On the possible parameters for a heavy Wigner matrix 

Proposition A.l. If a sequence {ak)k^i of real numbers is a parameter of a heavy Wigner matrix, 
then (afe_i)fc^i is the sequence of even moments of a Borel measure m with finite moments, i.e. 
for any k ^ 2, ak = J t'^''^^dm{t). In particular, if the parameter {ak)k^i is non trivial then one 
has flfe > for any k ^ 2. 

By the Hamburger's theorem |13] . a sequence of real numbers (^fi{k)^ is a sequence of moments 
if and only if, for any sequence {xk)k^Q of complex numbers with finite support, one has 

n{j + k)xjXk ^ 0. (A.l) 
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Let Xpf be a heavy Wigner matrix of parameter (afc)/j^i et let /i^v be the common law the sub- 
diagonal entries of \/NXn. As we do not change the parameter of a heavy Wigner matrices when 
we change sub-diagonal entries Xi,j into e|Xij |, where e is a random uniform sign, we can always 
assume that /ijy is real and symmetric for the task of the Proposition. Denote by (/^'•^•' (fc))fc^o its 
sequence of moments. For any sequence {yk)k^i of complex numbers with finite support such that 



yo = 0, we apply (A.ll with {xk)k^i = {N ^yk)k^i- we get 



= N ^ ai±kyjyk + o{l) 



where we have set aj, = whenever k is odd. Then, the sequence (afo^j^)fc^i satisfies (A.ll, which 
gives the proposition. 
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